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Abstract 

Techniques are presented for calculating directly the scalar functional 
determinant on the Euclidean d-ball. General formulae are given for 
Dirichlet and Robin boundary conditions. The method involves a large 
mass asymptotic limit which is carried out in detail for d = 2 and 
d = 4 incidentally producing some specific summations and identities. 
Extensive use is made of the Watson-Kober summation formula. 
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1. Introduction 

For specially symmetric situations it is often possible to relate the eigenvalue 
problem on a bounded region, Ai, to one on a covering space, the method of images 
being an example. As Pockels [1] has pointed out, in the case of the boundary 
condition, (n.V(^ — hcp^ Iqm ~ ^' ^^^^ technique is no longer applicable and one has 
to proceed directly. In classical physics such conditions arise in heat conduction [2]. 

In the present work we wish to analyse some cases of these Robin boundary 
conditions, evaluating the coefficients in the asymptotic expansion of the heat-kernel 
and also the corresponding functional determinant. The reason is that it is often 
useful to have explicit results for special cases in order to put restrictions on the 
general situation. In addition, Robin conditions enter into the mixed boundary 
conditions required for higher dimensional vector bundles. 

The domain of concern here is the Euclidean cZ-ball. The determinants on such 
regions have already been found, in the massless case, by a conformal transformation 
method but this technique is limited, at the moment, to d < 4 by the availability 
of the necessary transformation law. A more direct method is therefore needed. 

The recent determination by Branson, Gilkey and Vassilevich [3] of the C5/2 
coefficient is not sufficient to allow the extension to cZ = 5 as it is restricted to totally 
geodesic boundaries. 

Other calculations involving relevant material that should be mentioned are 
those by D'Eath and Esposito [4,5], Schleich [6] and Louko [7]. A parallel calculation 
by Bordag et al has recently appeared, [8] . 



2. The interval 



We work on / = [0, 1] and consider the simplest equation 



A0 



with 



dx 



■^^'^^ Tx 



This eigenvalue problem is discussed in Carslaw and Jaeger. We write the 
eigenfunctions 

(f){x) = N sin(Q;a; + 5) 



1 



where \ = and tan 5 = ajh^. The eigenvalue equation is derived to be 

tana + 
a — h^YiQ 



(1) 



The normahsation, N is easily worked out. 
If /j.Q = /ii = /i, (1) factorises into 

2h cos{a/2) + a sin(Q!/2)^ (^ ^"'^^^^^ ~ ^ cos(q!/2)^ = 0. (2) 

An eigenvalue equation similar to (1), but with the right-hand side reversed in 
sign, occurs in the theory of the three-ball with Neumann conditions at the surface 
(Rayleigh [9] Lp265, eqn.(2)). 

The eigenvalues are non-degenerate and, although easily computed, are not 
needed here since the Mittag-Leffler theorem can be applied 

(/iq/ii - z'^f-^ + (/ii + ho) COS z = (/iq/ii + K + K)W^['^- ^] (3) 
z \ (y, J 

a ^ ' 

where the product is over the positive roots. 

Incidentally, logarithmic differentiation, expansion in positive powers of z and 
the equating of coefficients allows one to evaluate the sums of inverse powers of the 
roots. (This technique goes back to Euler, c/ Rayleigh [9] II p. 279. Rayleigh does 
not differentiate but just expands the logarithms. See also Watson [10] 15.5, 15.51.) 

It is found that, for example. 



\_ _ /iq/?.! +3(/io + +6 

^P- QiK^i + K + ^i) 

1 hlhl + 6(/ig/ii + hohj) + 15{hl + hj) + SOhoh^ + 60{ho + h^) + 90 



90(/io/ii + /io + KY 
which can be confirmed numerically. Defining the (^-function by 



we see that the values Cl^)? ^ can be determined. 

One of our aims is to find the coefficients in the heat-kernel expansion. These 
are related to the residues at the poles of C,{s) and to the values C(~''^); ^ ^ '^^ • The 
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method we use here involves the time-honoured one of the asymptotic behaviour of 
the resolvent 



A + 

A 

as ^ oo. (c/Dikii,[ll], and [12,13]). 
More generally we can define 

A 

and it is straightforward to derive the large, positive m asymptotic expansions (c/ 
Dowker and Critchley [14], Voros [12], Elizalde [15]), 

r(.rr>^\ ^ \^ r^2xd/2-n-s r(g-c//2 + n) 

' n=0,l/2,l,... ^ ^ 



and 



^ fi;^ ^ (^-^^'^ " l/2)m-2"+i - Inm^) , 

^ n=0,l,3/2,... ^ ^ 





^1 + 




2m ' 






^1 


1 + 


2m ' 



and the coefficients C„ can be read off as 

^0 = 1, 



(7) 



If cZ = 1 it is sufficient to set s = 1 in (6) to get 

7./ 2x 1 V- r(l/2 + n)^ 

^ ^ n=0,l/2,l,... ^ ^ 

Comparison with (3) shows that z = im and that 

9n <^ , //, , o.sinhm /, , s , \ /„n 
R{m^) = -— In I (/ig/ii + m^) h (h-^ + Hq) coshm ) . (9) 

The large positive m limit is easy to find, 

~ — 1 + — + 

m + hi m + Iiq m 

1 

m 



Ci/2 - Vtt, (10) 
r(n + 1/2) 
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These check with known forms and we are now in possession of the values Cl"'^)- 
Our next objective is the functional determinant, 



D{-m'^) = exp ( - C'(0, m^)) (11) 



in particular -D(O). 
From (3) 



C'(0, 0) - C'(0, m^) - lnF(m2) - In (/iq/ii + K + h,), (12) 
where F is the factor in (9). We can thus write 

('(0, 0) = Jlrn^ ^C'(0,m^) + In F(m^) - In {hoh^ + /iq + /^i)^ • 

From (9) 

lnF(m2) ^m + lnm-ln2 + lnf (l + ^)(l + ^)\ (13) 

and using the asymptotic form (7) one finds that 

('(0,0) = -ln2-ln(/io/ii + + /io) (14) 

and so D{0) = 2{hQh^ + + ho). The same result emerges from the method of 
Barvinsky et al, [16]. 

As a check, expanding the logarithm in (13) gives the previous expressions, 
(10), for the coefficients. 

The Neumann (and Dirichlet) value for ('(0,0) is —In 2. But the Neumann 
case, = Hq = 0, has to be treated afresh since there are zero modes. The 
vanishing of the determinant, i.e. /^o^i + hi + Hq = 0, is related to the appearance 
of a zero eigenvalue but the corresponding mode does not exist. 

Turning equation (12) around, one gets the massive functional determinant, 

D(-rr?) = 2l(hQhi+m?)—^ h (/ii + /iq) coshm J . (15) 

Other, miscellaneous boundary conditions, [2] , can be treated in the same way 
with similar results. 
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3. Heat kernel coefficients on Euclidean balls. 



The higher dimensional interval, i.e. the ball, is treated by the same approach. 
The operator is the Laplacian and the standard mode expressions, which will not 
be repeated here, arc outlined by Moss [13]. (See also [17,18].) The eigenvalues 
A = Q!^ in (5) are given by the roots of 

J J a) = 0, Dirichlet 

; (16) 
PJp{a)+aJp{a) = 0, Robin, 

where P = h + 1 — d/2. The degeneracy for a given Bessel order is, for even d 

N^'^ = y^P\p' - 1) • • • (P^ - {d/2 - 2)% (17) 



while for odd cZ, 



The Mittag-Leffler theorem can again be employed. Since there are no normal- 
isable zero modes, the zero roots of (16) must be removed. We therefore write 



---^ = 711(1-^) 



(19) 



where ^ 

= ^4^' ^^^^^ {M^M- 

If /9 = 0, then p = is a special case, Fq — zJq — —zJi, and 

--^^ow=-in(i-5)- (20) 

Moss [13] was concerned to evaluate the heat-kernel coefficients, and we present 
here an equivalent treatment. The coefficients can be read off from the asymptotic 
expansion (6) with s = {d + 2)/2 ii d is even and s = {d + l)/2 ii d is odd. The 
requisite expression for the left-hand side can be obtained by repeated differentiation 
of the resolvent (4) or, equivalently, of the product (19), which is rewritten with 
z = im as 



m 



'w=7n(i+^) (21) 
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where is now given by 

FJm) — I Am), Dirichlet 

(22) 

= pip{m) + mlp{m), Robin, 

in terms of the modified Bessel function. (See also D'Eath and Esposito [4,5].) 
The required quantity is then (c/ [13]) 

= - ( - ' E 1- (m-^Fpim)). (23) 

The development of the logarithm follows from Olver's asymptotic expansion 
of the Bessel functions. In particular 

1 1 / °° U (t)\ 

In (m-P/p(m)) - -- ln(27r) + e - pin (p + e) - - Ine + In 1 + ^ (24) 

and 

In (m-P+^/^(m)) - ln(27r) + e -pin (p + e) + ^ Ine + In ("l + ^ ^^^"j (25) 

where e = ^/{p'^ + m^) and t = p/e. It is sometimes more convenient to use Olver's 
variable, z = m/e, and we will use the same kernel symbols for the polynomials 
whether they are written as functions of t or of 2:. 

The polynomials U^{t) satisfy the recursion relation [19] 

Un+i = -t')^ + ^ J\l - 5t') U^dt, (26) 

and the are related to the by 

Vn = U^-lt{l-t^){U^_, + 2t^). 



Equation (24) is relevant for the Dirichlet case. The Robin expansion is 

(27) 

V —1 P' / 

where ^ + l3tU^_i with Uq = 1. 



1 1 /■ 00 \ 

In {m-PFp) ^ e - pin (p + e) - - ln(27r) + - hie + In M + J] ^ j ^ 

^ n—l ^ ' 
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From the recursion (26) it is clear that U^{t), V^{t) and W^{t) each have a 
factor of f^. Extracting this factor, we write e.g. 

Unit) _ Unit) 
pn ' 

so that W^ = V^ + 

The systematic treatment of the summations in (24) and (27) involves the 
cumulant, or cluster, expansion, written generally as 

^ n=l ^ n=l 

with 

Sn = J2(-^)'~"(^-^y-Il% (28) 

where 5 = Yl^^i the summation is over sets of positive integers Ij including 

zero such that Yl^=i J — ''^ • 
For example. 



^1 


= Qi 


-Iqi 




S2 






S3 


= Q3 


— Q2Q1 






= Q4 


— Q^Qi 






= Q5 


— Q4^Qi 


— Q^Q2 + Q3Q1 



the number of terms increasing rapidly with the order n. 

Specific evaluation shows agreement with ref. [13], e.g. in the Dirichlet case, 

(24) 

= -r' + 

T,(t) = ^ + - HH,. + (29) 

^ 360 240 192 1152 ^ ^ 

^ 4 16 128 

^ , , 1 125 493 4 1621 ^ 82825 o 82825 ^ 

TJt) = + z^ z^ + z^ z^ + z^° 

' 1260 504 72 48 1536 3072 
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while, for Robin conditions, the corresponding polynomials are 

R2{M = Af + ^(2/3 + \)z^ ^ (30) 
^ X 1 l^s 59 o 1^2 2 259 4 7^ 4 1463 . 

Related expansions for the spin-half case were computed by D'Eath and Es- 
posito [4,5]. 

The next step is the differentiation in (23) followed by the sum over p, which 
is the hardest section of the analysis. 

To check the method, and for illustrative purposes, we firstly look at the Dirich- 
let 2-disc, q = 2 and Np — 2 — dpQ in (23). Algebra produces 

for p > 0. 

Since the bracketed term in (31) converges overall as p — > 00 it is legitimate to 
regularise it. Moss writes 

^-%lUlim^^-^ ' ' 



w? 2e / s^o V e^^m^ 26^^+1 
which introduces the series 



00 

C^(«.'»^) = E(;j^ (32) 



p=0 



and 



so that 



Then 



C,(.,m^) = X:^^^^^ (33) 
E (p2 = 2Ci(«. ^ ;;i2l ^ 2C,(.. m^). (34) 



(^+^) = ;! I2i (c.(«.'»VC.(.-i,™^)+^C,(.4. ™^)) . (35) 



Although (32) and (33) are very standard we proceed from first principles and 
write _ 

C.(.,m^) = J^y^ dtt^-'Y.Pe-^'^^-'^' 

which is related to the ^-function on a two-sphere. The asymptotic expansion of 
the heat-kernel, related to MulhoUand's result, is known to be 

Y.P^-'^'-^t E ^nf- (36) 
p=0 n=0,l... 

where 

Co = 27r, C„ = 27r^^S2„, n>0 

in terms of Bernoulli numbers. Of course the coefficients in (36) are best deduced 
by applying a summation formula to the ^-function, but it is convenient here to 
consider them as given [e.g. [20]). 

The asymptotic expansion of C2(s,m^) is then given by (6), 

,2-2s 



n=l 

Is 

+ ... 



f^(^'")~2(73T)+E("'^)' 2r(.) n! 



2s ^ 2s 



and, similarly. 



2(s-l) 12 



C.(.,™^)~f ^^^r.— . (37) 



In (35) the first term in brackets is immediately finite, 

C2(0,m^) = ----, 

and the other two can be combined, asymptotically, to a finite expression using 

(37), 



Ci(s + l/2,m2)-A^ (5_l/2,rTi2) 



rw r(,-i) 



2 v2r(s + i/2) r(s-i/2) 

0r 2s «r(s) 

lib 



2 (l-s)r(s + l/2) 

which equals 1/2 at s = 0. 



We therefore find, 

2 



and, further. 



Consider now the polynomial part of (31). The n = 1, 2 and 3 contributions 

are 

d \ Ti 25m2 11 



dm2 / e 32e7 IGe^ 
d 15m4 21m2 _^ 13 



rfm2 y e2 4^10 4^8 8^6 ' 
d Vt. 12155m6 11271m4 771m^ 57 



dm'^ J e3 512ei3 256eii 32e9 IGe^' 

After multiplying by Np , the summations can be effected using (34) and (37). 
Going back to (6), and writing out the first few terms, 

elementary arithmetic confirms the heat-kernel expansion derived by Stewartson 
and Waechter [21], and extended recently by Berry and Howls [22] to many more 
terms. We note that Stewartson and Waecliter's method also uses knowledge of just 
the asymptotic behaviour of the eigenfunctions, but the information is organised 
differently. 

The Robin case is similar to the Dirichlet one. The polynomials are, 
d _ 13 + 12/3 35m2 



rfm2 J e 16e5 326^ 

d Vi?2 8/32 + 16/3 + 15 24/5 + 9 21m4 



+ 3m' 



dm'' J £2 8e6 4eio 

d Vr. 20/3^+40/32 + 68/3 + 63 140/3^ + 532/3 + 917 
' — m 



dm'' J e3 16e7 32e9 

,3528/3+ 14217 « 16093 
+ "^^56?^ "^5T2?^- 

In two dimensions (3 — h. 
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The Robin coefficients, derived from a comparison with (40), agree with those 
found by Moss [13] and so will not be given here. There is no difficulty in carrying 
out the evaluation to any desired order. 

Waechter used the method of [21] to discuss the case d = 3, i.e. the three-ball. 
His result was corrected by Kennedy [23] who extended the calculation up to (i = 6. 

The present method also extends to higher dimensions without the need for 
any special summation techniques. Since the complications are simply arithmetic 
ones, symbolic manipulation can be used to advantage. 

A method has also recently been developed by Bordag, Elizalde and Kirsten 
[24] that allows any number of coefficients to be calculated in higher dimensions. 
Their method is similar to that of Barvinsky et al [16,25]. 



4. The functional determinant. 



In the general case, Voros [12] develops the connection between the two quan- 
tities D{z^) and A (2;^), where D is given in (11) and the Weierstrassian product 
is 



^(-')=n(i-?)-pi;r(4)', (41) 



A fc=l 



which amounts to subtracting the first [//] terms of the Taylor expansion of 
ln(l - z'^/X). Formally 



A 



(42) 



(See also Quine, Heydari and Song [26].) 
Voros' relation is 

<'(0, 0) = In - E FPCik, 0) V - E C«/2-. (1 + ■ • ■ + ^) V 

^ ^ k=l k=2 

where /j, = d/2 is the order of the sequence of eigenvalues and FP indicates use of 
the finite part prescription. When fj, = 1/2, (43) is equivalent to (12). 

From this equation it is clear that to find C'(0,0) it is only necessary to pick 
out the m-independent part of the right-hand side which can be done in the infinite 
m limit. We note from (6) that lnD(— m^) = — C'(0, m'^) contains no constant part 
in this limit and that the two summations in (43) are also not needed. Hence 

C'(0,0)= lim lnA(-m2), (44) 
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and this is our starting formula. Once C'(0,0) has been found, C'{0,m^) can be 
computed numerically from (43). A power series expansion in is easily derived. 



5. Determinant on the disc. 

As an example of the use of this technique consider Dirichlct conditions on the 
disc, for which /j, = 1. In this case (41) reads (the degeneracy is implied), 

A(-„.') = n(l + ^)-P(-^)=n (Pf^-m-'l,) e.p ( - jj^) , (45) 



where we have used Rayleigh's explicit results on the sums of inverse powers of the 
zeros of Bessel functions rather than the formal subtraction (42) . 
From (45) and (24), asymptotically. 



OO y 

InA(-m^) ~ ^Ar^^Mpln2 + lnp! -ln^7r + e-pln(p + e) 
p=0 ^ 

-lnv^ + J] 



Tjt) 



n 



4(1 +p)^ 

(46) 

We would like to do the sum over p and take the large limit. The integral 
representation of logp! is employed to give 



InA(-m^) ~-lny(27r) + m- ^lnm+ A''^^^ 

p=0,n=2 



12^ P 2e3^ ^ 

p=0 p=l 



phi — - — \- e — p (47) 



Ap 2 p \2 t e^-l) t 



dt 



where some initial terms have been extracted for convenience of doing the summa- 
tion and taking the asymptotic limit. It is readily confirmed by expanding in m/p 
that the sum over p converges. 
Consider firstly, 

lim > — 5— 

m-^oo 24 e'^ 

which can be evaluated using (37) to give 5/12. 
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Next look at 



^ I fl 1 1 \e-*P , 1 



p=l 



The sum over the integral diverges because of the behaviour at t = 0. Adding and 
subtracting the offending term, and performing the sum, yields 

Jo U t 12^e^)tiet-l) ^T2^^p 
Dealing with the final summation first, we make use of the standard formula 

[27] 

J2 (---)= + ^ - 2 ^ Ko(27rmp) + 7 - log 2 (48) 

p=i P ^ ^ p=i 

the last two terms of which are the required constant contribution. The integral 
term can be regularised as 

The remaining terms in the summation in (47) are now evaluated in piecemeal 
fashion. 

The ln(e/p) term in (47) can be found by differentiating the particular Watson- 
Kober formula [28,27] 



oo 



1 1 OS (1— 2s)/2 °° 

(ni2s) = E (-L-i-) + ^n(^-^)/^i^(,,_,),,(2.nm) 

f=i (50) 

mi-2*r(s-l/2)0F 

2~ ^ 2r(s) 
valid for Res > —1/2, in particular near s — 0. We find 

oo 

^(Ine-lnp) = — - In ^27r) - - Inm + - In (l - e'^'^'^) . (51) 

Alternatively, the asymptotic part only of this follows by differentiating 
°° 1 1 - 1 

E(^-^) = Ci^(2^)-Ci(^) + ^ 

13 



and using (37). The same goes for (48). 

To treat the (e — p) term, twice the following term, i.e. m?/2p^, could be 
replaced by m^/2e^, plus a known correction, using (48), so allowing the earlier 
result (35) to be applied. However it is preferred here to use directly Watson's 
general formula. 



1 ^1 rn^n N _ ^Y{s - l/2)mi-2- 1 



^ ( (p2 + ^ V ^ y p2s+2n / 2T(s) 2 



- -m 



+ fS™^'"''^^' 5:p(^^-^)/^K(,,_,)/,(27rmp) (52) 

M-l 



E 'W^-C^(2. + 2n) 

n=0 V ^ / 



which is valid for Re s > 1/2 — M and plays an important role in the sequel. If M 
is set equal to 1, one gets (50) and the s — > 1/2 limit then yields (48). We now need 
to take M = 2. This gives 

1 , \_ 0Fr(.-i/2)m^-^- 1 

^\e^' +^^2^+2^ 2r(s) 2 

+ f;p(--)/^K(,,_,)/,(2.mp) ^''^ 



,2, 



-Cii(2s) + WCji(2s + 2) 
valid for Res > —3/2. In particular at s = — 1/2 

(e-p ) = — (21n2 + l)--7m2--m--m^lnm^ > — 

^ ^ It,' ^2' 2 4 TT^ » 12 

(54) 

which, in conjunction with (48), provides an alternative derivation of (35) and (38). 
Incidentally, the integral representation 

/CXD 
e-^\t^ -Xy-^l'^dt (55) 

allows the summation to be performed, 

_my. Ki(27rmp) T \ {t^-xyi-^dt 



p=l 

m 

TT 



/•oo . 

A ^(^2- 1)1/2 



(56) 



dt. 
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Equation (54) then agrees with a formula derived by Ehzalde, Leseduarte and 
Zerbini [29] in a different connection by a longer method. 

Harking back to (47), all terms have now been dealt with except for 

^ V p + e 4p J 



If this expression is differentiated, (48) and (54) employed and the result inte- 
grated we find, 

Efpln-^^ + = -(m^ + -) Inm^ + \m'^(-/ - In 2) 

V p + e 4p J 8^ 3^ 4 ^' ' 

(57) 

m^Ki(27rmp) 1 ^ K^{2T^mp) 1 h ( 

+ 7^ p +2^^ p^ Y2^''^"^«^"^^- 



p p 



where the constant of integration has been found by setting m to zero. This is 
allowed because (57) is exact. From (54) 

lim^y^l(^ = l (58) 

p=l 

which can be derived from (56) and also by replacing Ki{z) by its leading term, 
~ as ^ — > 0. 

Replacing Kq{z) by its leading term as 2; — > 0, i.e. — ln(2;/2) — 7, leads to 

p=i (59) 
= C^(-l) + ^(ln2-l-lnm) 

where we have used 

Sc;?(2)=Cfl(-l) + :;^(ln27r + 7-l). 



27r2^"' ' ' 12 

These results enable equation (57) to be completed. 

As an aside, a useful expression arises from the Fourier transform representation 

K^z) = ^-^l\2zlarV{y + 1/2) ^ ,2y+i/2 ^t- (60) 
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Thus ^ 

KQ{2Tvmp) -sr^ cos pt dt 



The sum over p gives a BernouUi polynomial in the range < t < 27r and 
splitting up the integral into 27r lengths we find, after rescaling, 

^ ~ i-i ^ {{t + + 4m^) V2 ^^^^ 



where 



52(^) = 



12 

The integrals can be performed, leading to the summation 

^ 24 ^£5j^3^ [ ^ V ((n - 1)2 + 4m2)^/' + n - l/ 

3((3n - 1) v^(n + if + 4m2) - (3n + 1) vT('i - + 4™^)) 

(62) 

which is convenient numerically as a means of checking the foregoing expansions. 

These results can be generalised to sums of the form Ylip^oi'^P) /P^'^ • We give 
only the formula 

P=l (63) 
(_l)n+l22n-1^2r ' 

" {2n)\ 
and also note that 



(^S2„(ln2 - lnm + 7 - V(2n)) + 2nCfl(l - 2n)^ 



lim — ^^'^^Trrnp) ^ ^^+;_2 (/ - 1) 

m^O TT" ^ pn (^ + Z) 



when / + n is even. 

Incidentally, taking the m — > limit of (62) and comparing with (59) gives an 
expression for Ci?(~l) 

oo 
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Numerically, from (65) we find Cfl(-l) ~ -0.165421143700045092139. A summa- 
tion for Cr{~^) is given by Elizalde and Romeo [30]. Similar summations can be 
found for C^(l -2n). 

All parts of (47) have now been evaluated. Collecting the constant terms we 
arrive at 

C2baii(0,0) = 2C^(-l) + lln2 + ^_lin7r«. 0.773714, (66) 

in agreement with Weisberger [31]. 

As a tactical point we note that the equations we have developed contain 
far more information than needed for the specific purpose of picking out the m- 
independent part and so, from this point of view, the method is not very economical. 
A more streamlined version is presented in section 8 where only the essentials are 
retained. 



6. The Robin case. 

In this case the appropriate expansion is (27) and we firstly need ^l/o:^. It 
is straightforward to derive, 

^^ = 4(^Tl)(^ + ^) ^g^^ 
= /3 = 0,p = 0. 

We remark that Lamb's formula, [32] 

1 _ 1 

is not enough to prove results like (67) but does provide a test when (3 = p. A 
rough, numerical check of (67), and similar formulae, can be obtained in the p = 
case by using the roots given in Table III of Carslaw and Jaeger [2]. 
Instead of (45) we have 

A -m^ = 11 F„) exp - — -(l + ) 
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and so 



InA(-m^) ~ ^ [pln2 + lnp! -ln-\/(27r) - ln(/? + p) + e - pln(p + e) 



Corresponding to (47) we find (/? 7^ 0) 



lnA(-m2) ~-lnv^27r)-ln/3 + m- -lnm+ ^ N'f^^ 



2 

p=0,n=2 



rr)^ 2 If 

which can be reorganised in the same way that (47) was. 

The only sum not covered by previously given expressions is 

which, using (48), can be converted to 



■ 2p 

p in \- e — p 

p + € 



(68) 



(69) 



Putting all the pieces together and selecting the m-independent part yields, if P ^ 0, 

CL(0, 0) = 2C^(-1) - ^ - ^ ln2 - i luTT + 21nr(l + /?) + 2/?ln2 - In/?. (70) 

The case when /3 is zero must be treated separately because of the appearance of 
extra zero modes. The p = term contributes an additional In 2, which is clear from 
the fact that Jq{z) — —zJi{z) ~ and the result is the Neumann determinant, 

^(0, 0) = 2C^(-1) - 1 + 1 ln2 - 1 luTT. (71) 
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7. The four-ball. 



We will pursue the calculation of the functional determinant in detail only for 



B^. In this case the order /j. is two and equation (41) reads 
A(-™=)=n(l+^)exp(-^ + ^) 

p,a P P P 



n (p!2^m-^/J exp ( - 



m 



+ 



m 



(72) 



4(1 +p) 32(l+j9)2(2+j3) 



) 



and (46) is 



InA(-m^) ~ ^p^(pln2 + lnp! - In v/27r + e - pln(p + e) 



-lnyi+ J] 



+ 



4(1 +p) 32(l+p)2(2+p)y 

(73) 

The situation now is that more terms in the summation over the polynomials 
T„ become relevant. 
(47) is 

Sm^ 5m'^ bmP w? 



p=l,n=4 p=l ^ 



+ 



+ 



1/1 1 



+ 



360 V e P 



24e 24e3 166^ 16e6 4e2 
ISlrn^ 1349m4 2431m*5 1105m* 



720e3 



■ 2p 

pin \- e — p 

p + e 



960e5 

,2 



+ 



11526^ 1152e9 

4 

m 

4(1 +p) ^ 32(l+p)2(2+p) 



(74) 



1 e 1.1 1. 
-In - H ( ) 

2 p 12 > 



+ 



1 



2 t 12 ^ 720 ' e* - 1 y t 



+ 



,-tp 



dt 



Rearranging the summations of the Rayleigh terms for future use, 



E 



m 



+ 



m 



4(1 +p) 32(l+p)2(2 + p) 



^ OO 



1 

OO r 



m [l — p h 



1 



1 



+ 



p p{l+p)J 8 \p p{l+p) p'^{l+p) 



-T 



m ( 1 — p — 



H m 

8p 



4 

4 



TT 



-4 



(75) 
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The convergent summations would be needed for a finite mass calculation but 
are irrelevant here and (75) can thus be replaced by 



4 ^ 



m \ p — 1 + 



P 



1 

8 V 



+ 0{m^). 



The formulae needed to analyse (74] are extensions of those used earlier. We 
set M equal to 3 in (52), to get 



^(l l_ 

2^ [ /Is „2s 



m 



2 



s(s + l) \ _ 0Fr(s - l/2)mi-2* 1 



- C^(2.) + sm\^{2s + 2) - ^^^m^CH(2« + 4) 

(76) 

and then let s tend to —3/2, yielding 



(e — p p )= — Cr(— 3)H — m m m mm m7 

^ 2 8 » ^ >>Kv / g 2 16 8 

+ -m^(41n2 + 3) + ^^jL_J^i. 

(77) 

Another necessary equation follows upon differentiating (53) with respect to s 
and then setting s = — 1, 



V (2e^(lnp - Ine) + m^) =m^ Inm^ - \-n? - In? In ^/(27r) - 

-72] ^5 ^E^a 2C^i(-2). 

(78) 

As a small check, we note that Ch(~2) = — Cij(3)/47r^ and the right hand side 
correctly vanishes when m is zero. 
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The final equation we need is 



—(37 -31n2 + l) + — m^ + —m^ 



32p y 32 ' ' '12 48 



3m ^ KjpmTTp) _ 3 ^ JCo(2m7rp) ^ 1 , ^ _ /_on 



which is obtained in a similar way to (57) and involves using (63) and (64). 
The integral is again regularised 

s^oJq \2 t 12 720 e*-iy dt^ - 1 ^ ^ 

The power terms are most easily dealt with by partial integration which gives 

tn'-'^-r—^dt={-iyT{fi + s)Cn{f^-j + s). (81) 

For the specific case here, j = 2 and the finite part of the power terms is 

a-3) + ici(-2)-^7+3^. (82) 
The remaining piece of (80) is 

Jo e*-l 



00 

4.S-1 



2 3 1 

(e* - 1)4 + (e* - 1)3 ^ (e* - 1)2 ' 



(83) 



These integrals are examples of a particular Barnes zeta function, 

, , , iT(l-s) f e^('^-")(-z)«-^ , 
C(s,«) = 7^ / — T\7 

00 



27r (e- - 1)' 



E/n + r — 1\ 1 
\ r — 1 y (a + n) 



(84) 



n=0 

In the present case a = r and the lower limit can be adjusted so that 



n=l ^ ^ 
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Expanding the binomial coefficient as a polynomial in n using Stirling numbers 

'n-r 

r-1 



I V S^''^ n!"-^ 



gives 

CM r) = ^r'^ -k + 1) (85) 

^ k=l 

in terms of the Riemann ^-function. 
Explicitly for d = 4 we need 

Csis, 3) = ICr{s - 2) - ^Cr{s - 1) + Cr{s) 

CAs, 4) = IcRis - 3) - CRis - 2) + ^CRis - 1) - CRis). 

Of course, one could legitimately leave the answer in terms of the Barnes 
^-function and the multiple F-function. 

The finite part of (83) is evaluated at s = to give 

2CK0, 4) + 3C^(0, 3) + (2(0, 2)-7(2C4(0, 4) + 3(3(0, 3) + (2(0, 2)) = 

^C^(-3) - ^C^(-2) + ^C^(-l) + ^7 ^^^^ 

and this has to be added to (82) to get the total contribution from the integral 
term. 

With these and the earlier equations it is possible to evaluate the summations 
over p in (74) and to extract the m-independent part. A point to note is that the 
convergent terms coming from the T3 polynomial all contribute to this constant 
term using (37). 

Regarding this, a minor calculational point should be remarked. Because of the 

factor, the sum over p can be extended to p = when convenient. This means 

(2) 

that a factor of Np ' /2 can be inserted without penalty so allowing the simple 

asymptotic form (37) to be used. A more refined approach is given in section 8. 

Combining all m-independent contributions by hand we obtain 

1 1 1 1 1 73 

C4baii(0) = 3C^^(-3) - -C'(-2) + -C^^(-l) + ^1^2+ ^ 0.002869, 

agreeing with a previous result [33] . 

We note that there is cancellation of the terms in (82) with those from the 
summations, leaving just those in (86). The 7 terms always cancel. These facts are 
generic as will now be shown 
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8. The d-Ball. 



With the experience gained in the previous sections, the calculation can be 
readily extended to arbitrary dimensions. 

The degeneracy is expanded as usual, (for d > 2) 



d-2 



which is an odd polynomial for d odd and even for d even. The coefficients are 
Stirling numbers. 

For even balls, it suffices to discuss 

OO y 

InA(-m^) ~ ^p"^" (pln2 + lnp\ - In V27r + e - pln(p + e) 

-ln^^+f;^-Ray(m,p)). 



OO 



(88) 



= V/" pln^ + e-p--ln- + 
^-^ I p + e 2 p ^-^ 



n=l 

dt 



where Ray(m,p) are the 'Raylcigh' terms. 

It is enough to decide what terms are irrelevant for the m-independent part of 
the large m limit and we first note that the Rayleigh terms are so irrelevant. Thus 
it is not necessary to find the sums of inverse powers of the eigenvalues, interesting 
though these may be. 

Again we go through the various contributions. The analysis is somewhat in- 
volved and fragmented, but not without interest. Much of the difficulty is bookeep- 
ing. 

First consider 



= (-l)"-^m2"-'^(e2^+" -p2«+«) + ^ V-l)^-W^- 

q=0 q=0 

(89) 

for a = 1. 
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From (52) the relevant part of J2p (e^^"*"^ -p^^"""^) is just -(^(-2^' - 1) and so 
that of ^pP^'^i^ — p) equals —C,ji{—2v — 1) since terms involving m can be ignored 
and the second sum in (89) is ultimately irrelevant. 

The next contribution is 

-Jy/Hn- 
2 p 

p 

whose relevant part, —('j^{—2h')/2, is obtained by differentiating (52) and setting 
s = —u, after using = p^ + m?. 
Now consider 

^p-+Mn-B^ (90) 

P 

which is the hardest sum to evaluate. Our method involves differentiation with 
respect to m^, use of the Watson-Kober formula (52) and then a back integration, 
the relevant term emerging as a constant of integration. 

It is clear from the examples in previous sections that only the Bessel functions 
in (52) ultimately contribute to the relevant term. However we write out (52) in 
general form, setting M = q + 2 and s = —{2q + l)/2 {q e Z) to give 



-Cj,{-2q-l) 



_ 1 _ rfa +l/2)(l + 2,) , _ ,„ 2 ^ g ( 1 1 )A 

2 20F (9+1)! V k' ) 

I J r(-9-i/2) ^ 

(91) 

only the last term of which is actually needed. 

Returning to the required contribution, (90), sufficient members (v in number) 
of the Taylor expansion are subtracted to render the summation finite. We denote 
this regularisation by Yl* . Then we find, using (89), 



p q=0 \ y / p 



2ui2 Z-^ V ^2 ^ v' 

p r- 

(92) 



where ^*(e^^+^ _p2q+i'j j^g^ ^j^g left-hand side of equation (91). 
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Integrating (92) we see that all terms in (91), apart possibly from the Bessel 
function ones, will produce irrelevant contributions and so the problem is reduced 
to evaluating the massless limit of the summations 



oo „ 
13=1"' 



pl+q ^ > 



where q ranges from —1 to v. 

The integrations can be performed using a recursion given by Watson [10] which 
reduces the 2v exponent to zero in steps of 2 finally allowing use of the integral 



/ 



The result of the recursion is given in terms of terminating hyperbolic Lommel 
functions [18], 7.14.1 (7), 

J x'^''-'iK^^q{x) = -2vxS2y-q-i^q{x)K^^q{x) - xS2y-q^q^i{x)K ^{x) . (94) 

Written out, the Lommel functions are 



^2u-q-l,q{^) — ^ ^ ^1" 



x"^ 

42(^_l)(^_2)(^.-g-l)(^.-g-2) 
^ X^ ^ 



and 



ev{y-l){v-q-l){v-q-2) 
x^ 

Since the massless limit has to be taken, some simplification can be made at 
this point. The leading term in Kq{x), g > 0, as x — > is 2^~^(g — ly./x'^, while 
Kq{x) ~ — (ln(x/2) + 7). For g = — 1 and = it is easily checked that the Kq 
terms in (94) go out as a; — > and that, for < g < (z/ — 1), so does the term. 
The explicit form of the Lommel function then yields 

lim /a;2^-«Ki,Ja;) = -22^-«-^(zy-g-l)!zy!, -l<q<u-l. (95) 
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When q = the inverse power terms in the series expansion of the Bessel 
functions cancel on the right-hand side of (94) leaving just the logarithm terms and 
positive powers which yield the limit 



^ fe=i ^ 



(96) 



To regain (93) we set x = 2Timp. As m ^ 0, (96) diverges as Inm but this will be 
cancelled by other, irrelevant terms that we do not need to calculate (except as a 
check) . 

Up to this log term, we find that the massless limit of the Bessel terms obtained 
by substituting (91) into (92) and integrating with respect to is 



^'(-^)'c.(2..2)x: ^-^r 



27r2-+2 ■'Kv ^^^^r(-5_i/2)(5 + l)!(5_^) 



((In. + 7 - ^ Ej) + 2) - C;.(2. + 2)) 



(97) 



which, in view of the identities. 



^ ^ (2iy+l)!r(-iy-l/2) 



and 

lu + iy. ^ u-q \ q ) / [ V ) ~ f^/c 
simplifies considerably to 



(21^ + 1)! v — q \ q ) \ v I ^ k ^ k' 



Cr{-2v - 1) + / ,, i?2.+2 In 2 (98) 



(-1)'^ 
2(zy+l)" 

with no fractional part. The ultimate simplicity of this result suggests the existence 
of an easier route. 

When integrating (92), the constant of integration necessary to give a zero 
value for the regulated sum, In (2p/(p + e)), in the massless limit is the 

negative of (98) and, finally, this is the relevant term coming from (90). 

The remaining contribution in (88) is 
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The first thing to note is that for n > 21/ + 1, the sum over p of p^^T^/e^ 
converges and the result is irrelevant. This can be shown by expanding the binomial 

= (e^ — m^) ^ and noting that, according to (37), only terms like vff je^^^ sum 
to a relevant part in the asymptotic limit. Alternatively, the limit, [13], 

r(r + l/2)r((n-l)/2) ^i nnn^ 

gl^T- 2r(r + n/2) ' ^ > 0' ^ > ^ » 

can be used. 

If the bracket in the integrand in (99) is completely expanded and the integra- 
tion performed, one gets the standard asymptotic formula, 

2 e'-l) t f^^^ ' 2k{2k - 1) j92fe-i 

f-rjl) 

n=l 

where the last equality comes from Olver's relation between the U^{1) and the 
coefficients in the Stirling approximation, [19] (2.15), (2.16). This shows that T'^(l) 
vanishes when n is even and that T2i^_i{l) = (— l)'^-B2fe/2/c(2A; — 1). 

If the n summation is restricted to the possibly relevant n-values and appro- 
priate terms added to and subtracted from the integral, one finds for the possibly 
relevant part of (99) 

'Tjt) r„(i) 



n=l 



pn 



(101) 

roo/-, .2k-l -1 \ p-tp 1 

G-^E<-)'-^w^?^)v4 

The two components of this expression are looked at in turn. Consider first 

Ep"E(^-^)- (m 

p=l n=l ^ ^ ^ 

The Tj^{t) are polynomials in t. Set 
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Hence one has for (102) 

oo 2iy+l , s oo 2t/+l ,2 ^2\v 

J: T„(l)(e= - mr(^ - ^) + E E ^^^^I-nW- (103) 

p=ln=l,3, ^ ^ p=l n=l 

Most of the terms in the second summation can be dismissed as irrelevant - if 
they are divergent because of the m-dependence in T' and if convergent because of 
(37). The only relevant contribution is when n = 2^ + 1, 

^^-+^^^^ (104) 

A convenient computational way of finding the asymptotic limit of this sum is 
to extract a factor of (1 — t^) — w? /e^ from T^, expand the resulting polynomial, 

and then use formula (100) to get 



y L_W(^ ^ y ^ 2.+i)^2 y t ^ y _JH 

p=l i=Q p=l i=0 



(105) 



Jo 

In the first summation in (103) we can write, using (89), 

f (.^-m^r(i--L) ~X:(^)(-ir-'"''''-''E*(^''-"-p''-")- w 

For nominally divergent terms to be relevant, there must be no factors of m. 
(Refer to (91).) This means that u = q and then all n values are divergent. These 
yield 

E T.iDE* ~ 

n=l,3,... 

v 

T2,+i(l) (In 2 - 7) - 5]T2fe_i(l) Cni^k -2iy-l), 

k=l 

(107) 
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since the relevant part of X^*(e^'^~"' — p^^""^) is —Cain — 2i/), for n < 2^ + 1, 
according to (91), and ln2 — 7 if n = 2^- + 1 according to (48). 

For n = 2^ + 1 there are also convergent terms in (106) which can be evaluated 
using (37) to give a relevant contribution, 



^-«(i)e(:)(-i)'E 

0=1 ^ ^ ^ p 



P 



The combination of (107) and (108) is 



(108) 



+ Wl)(l-2-7+iE^ 

^ Q=l ^ 



-1)1^1^1 



g)!r(g + l/2) 



(109) 



The integral in (101) remains to be evaluated. Performing the sum over p and 
introducing a regularisation, it becomes 

lim / --- + y(-l)^Bofc^— + — ^ t^-^ , / dt. (110) 

The power terms have been dealt with earlier in (81). The remaining part is 



Writing 



with 



poo -, 1 

lim / -A-rt^-^ " , -, dt. (Ill) 



^ ' dP e* - 1 ^ (e* - 1)' 

= Mp-') + (/ - l)4i-'\ 2<l<j 

A^l,=j\ (112) 

0') 



Ai^^ = 1 

(which is easily iterated by machine) we find for (111) 

,(2i/) 



Ji^o^"^' ^^"^^'+1^"'^ + ^^ ^^^^^ 



/=1 
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the finite part of which is 

E A^'^^ (C^+i(0' ^ + 1) - 70+i(0, l + l)). (114) 

Equation (85) allows the Barnes ^-functions to be rewritten in terms of the 
Riemann (^-function. Combining the two parts of the integral term gives 

lim (It{sKr{s - 2u) - r{s - l)Cfi(s -2^.-1) 



- E(-l)'^'7fSir(. + 2k- 1)Cr{s + 2k-2u-l) 



k=l 



2u+l l+l q(k) s 

+rw j:A''""%icH(--*^+i))- 

1=1 k=l ' ^ 

The cancellation of the individual divergences is a check of the analysis and 
implies the identity 

^^^-'^ (2! + 2)1' ^^'" + 1) + E E A^"^ Y U-k + 1) = 0, (116) 

where k will actually be even. 
The finite part is 

ICni - 2u) + CRi-2u - 1) - (7 - l)U-2u - 1) - l^-^f^n^i^ + 1) 
- E(-l)'^'|iytr(2^ - l)Ci?(2fc -2u-l)+ T2,+i(l)^(2z. + 1) 

2u+l l+l q(k) 

+ E E 4'^' ^ i(n(-k + 1) - 7C^(-'^ + 1)) 
/=i fc=i 

which can be reduced using (116) leaving, 

ICr{-2i^) + Cr{-2i^ - 1) + CR{-2iy - 1) + T,^+iil) {^{2y + 1) + 27) 

- E(-l)'^'7|Yi^(2/c - 1) Ui2k - 2. - 1) + E E ^f"^ Y + 

(117) 

A certain amount of simplification occurs if (117) is combined with (109). The 
first summation in (117) goes out and a 7 term cancels so that the i\) can combine 
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with the last 7 term to remove any remaining 7-dependence. Including also the 
relevant parts of YlpP^'^i^ ~ P) — IX^pP^'^ ln(e/p) obtained earlier, we find 

2u+l l+l q(k) 

1=1 fei ■ jjjgj 



Added to (104) and (98), this would be the final answer for (88), but first we 
manipulate (118) a little. Specifically, the order of summation is changed, 

1^1 k=l 



= E C'^{-k) + ^C^(-2- - 1) = E C'ni-k) 

k=l k=l 



2u + l 

where 



Mr - E ^r % 

l=k 



(r) 

We have used the Stirling number value Sr = 1- We have also used the values 
Si^^ = (—1)'/!, and the sum rule. 



1=1 



(which is easily proved from the recursion relation (112)) to alter a lower summation 
limit from /c = to A; = 1 in the intermediate algebra. 
(118) now becomes 



^{2v) 



C^,(-2^-l)+E<'^^C^^(-/^) 
fe=i 

+ Wl)(ln2 + El + i:^ 

^ k=\ q=\ 



, (119) 
g)!r(g + l/2) 



The final contribution to (88) is the asymptotic form (105). Adding this to 
(119) and subtracting the Bessel function contribution (98), we obtain the complete 
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Dirichlet expression, 



E ^^"^ c^(-^) - ^ ZlV 2 + /' {t) dt 

k=l ^""^^ -^0 



(120) 



The particular values 



are readily obtained and give 

Also M^l% = V(2z^ + !)• 

The final step to finding C'(0) on an even ball is to compound (120) with the 
expansion of the degeneracy (17), 

^p'^ = Td^.i^^^^^P"^' r=l + (d/2-2)^ (121) 

The explicit evaluation for any specific dimension d is straightforward, all but 
the rational part being quickly done. The integral term in (120) is the hardest to 
find since it involves the cumulant expansion, (28). 

The examples that have been calculated agree with those exhibited by Bordag 
et al [8] and in [33] for (i = 4. The d — Q expression is repeated here for completeness, 

4027 1 1 
^6baii(0) = -6486480 " 756 ^""^ + 60^^^"^^ 

- ^Cfl(-4) + ^Ck(-2) - ^C^(-l) 

fti -0.000392. 

We also give the coefficients of the Ci?(~^) terms for d = 8 in the form of a 
vector labelled by A;, 

1 7 79 5 31 1 19 



135' 144' 2160' 72' 648' 30' 990 
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9. Robin conditions. 

The Robin case proceeds in very similar fashion. The starting point is the 
expression 



lnA(-m2)~f;p2^Ln^ + 6-p+ilni-ln(l + /3/p) + 5^i^ 

p=l n=l 



- Ray(/3,m,p) + 



11 1 \ e-*P 



2 t -1 J t 



(122) 

Again, so far as relevant terms go, the Rayleigh contribution can be ignored, 
as can the n > 21/ + 1 terms in the summation. 

Prom the formulae in section 3, it is not difficult to derive the relation 



■n=0 

which implies (c/ [8]) 



i?„(/3,l) = Tjl) + 



1) + 



n 



as 



This allows one to add and subtract terms to rewrite the effective part of (122) 

lnA( - m^) y^V^'lpln^^ + e - p+ - In - - ln(l + /3/p) ] 
V > z-^p y p _|_ e 2 p / 

+ E f f E ^-t'^-*);^"''^'') + fl„(/3, 1) (i - 

p=l L n=l n=l ^ ^ ^ 



1 1 



t 



2k-l 



+ 



k=l 



{2k)\ - ly t 



-tp 



dt 



(123) 



The only really new term is ln(l + (3/p) which is easily reduced to 



^V'^ln(l + /3/p) = 



7 _^2.+l 



/3 



/32'^V(1 + /5) (124) 



2j^ + 1 

Most of the other terms are identical to the Dirichlet case but the cancellations 
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are not as complete and we find for the total expression corresponding to (120), 



k=l 



+ 



1 1 

k 4u 







Jo Jo 

^^2/c-l^^^^^ ^^ + ^^2^+^^^^+2^.+ 1^^2g(^.-g)!r(g + l/2)- 

(125) 

Again this has to be compounded with the expansion of the degeneracy. 
For d = 6 we find 



CM = 



9479 



32432400 315 



11 ^ 517 ^2 41 ^3 



15120 



_ 1^5 

1512 160 45 
1 



7^ + T^^'^^ " ^^'^] + 1^ / " ^^^^^ + 

+ ^C^(-5) + ^C;^(-4) - ^C^(-2) - ^C;^(-l) 

again agreeing with [8] after some rearrangement. 
The coefficients of the Cfl(~^) for d = 8 are 
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31 1 19 



135 ' 144 ' 2160 ' 72 ' 648 ' 30 ' 990 



10. The massive case. 

Although our main interest is with massless fields, we indicate in this section 
a means of computing the finite mass determinant, lnD(— m^) = — C'(0, m^). This 
will follow from the basic formula (43) which is repeated here 

C'(0,m2) =C'(0,0)-lnA(-m2) 

~2fc M 1 ^2fc (12Q) 

+ FPC(M)V-E^./.-.(i + - + itiT)V 

lnA(— m^) can be found from the properties of the Bessel functions but the finite 
part 

FPC{k, 0) = C(A;, 0) (A; not a pole) 

/ R\ (127) 
= lini f ({k + e, 0) j (A; is a pole). 
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is a transcendental quantity and its direct computation would involve an analytic 
continuation of the ^-function. However, since we know, in principle, the precise 
large-m dependence of the other quantities in (126), FP({k, 0) can simply be read 
off. This procedure is now illustrated for the disc, [//] = 1. 

From the known asymptotic behaviour of C'(0, w^) in terms of the heat-kernel 
coefficients, and from the dependence of the quantities in (47) obtained in section 
5, one can confirm that the only term in In (A(— m^)/Z)(— m^)) that does not cancel 
is proportional to m^. Therefore, since ({s, 0) has a pole at s = 1 of residue 1/4, 

li2l(«''.o)-4(J^)=|(^-i-'"2)- <i2») 

On the disc, inserting the power series for 1^ into (45), 
In A(-m2) = V Ar(2) \\n(l + V ^ - 



(129) 



with the extension to higher dimensions being obvious. 

The Rayleigh terms can be found from the Taylor series (cumulant) expansion 
of the logarithm and, making this expansion, one obtains a power series in vn? 
that begins with the first term not subtracted. Then, in terms of the well-known 
Rayleigh functions defined by 

(129) becomes 

OO OO 

1=2 p=0 1=2 

In d-dimensions, the Z-sum starts at [d/2] + 1. 

The expressions for the Rayleigh functions allow the sum over p to be per- 
formed. The first few terms are 

C'(0, m') = C'(0, 0) - ^(7 - 1 - In2)m2 - ^{27t^ - 15)m^ 

-^(35-67r^ + 24C^(3))m« + ... (^^O)) 
^ C'(0, 0) + 0.557966m2 - 0.0246834m^ + 0.0020103m^ + . . . 
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Of course, this series follows more directly by expanding the definition (5) but it 
would then be necessary to find the remainder (128) by an independent calculation. 
For comparison, the large positive m expansion is 

1 TT 1 TT 2 STtT 

C'(0, ^ -™^(ln™^ - 1) + -™ - - In™^ + — + 3^ + ^ + . . . (131) 

and a graphical analysis shows very good agreement between this and (130) for 
in the range 0.6 to 0.8. We should remark that Elizalde [15] has performed 
numerical calculations on similar expansions. 

11. Comments 

In this paper, which is a mainly technical one, we have presented methods for 
the evaluation of the functional determinants on the d-ball for Dirichlet and Robin 
boundary conditions applied to a massless scalar field. 

The information required is the asymptotic behaviour of the eigenfunctions, the 
eigenvalue properties being introduced via the Mittag-Leffler theorem. The method 
involves extracting the mass-independent terms from a large- mass asymptotic limit. 

For the disc and 4-ball, the calculation is done in more detail than required. 
In these cases, complete expansions are obtained which allow the necessary cancel- 
lations of the mass-dependent terms to be confirmed. The particular summations 
encountered in these explicit evaluations may have independent interest. 

For arbitrary dimension (actually only even here) a systematic procedure has 
been presented and an explicit formula given for the functional determinant in terms 
of coefficients evaluated by recursion. 

The method of [24] has been also used by Bordag, Geyer, Kirsten and Elizalde 
[8] recently to evaluate the functional determinants on balls directly from properties 
of Bessel functions in an effective and systematic fashion. There are many similari- 
ties in the intermediate expressions to our work but the details are different. They 
give precise results as far as d = 6 and indicate a general method. 

The spin- 1/2 values are available for d — 2,3 and 4 (J.S.Apps) by conformal 
methods and one anticipates only technical differences in extending the present 
approach to this case. 

The corresponding calculation on spherical caps should also be possible. (See 
Moss [13] and Barvinsky et al [16].) 
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